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Abstract 

A new mechanism of double-component convection is discovered. It emerges in a 
horizontal layer of Boussinesq fluid as a stable stratification due to flux boundary 
conditions is added to an unstable gradient specified by fixed boundary values. 
Driven by this mechanism, steady finite-amplitude flows reminiscent of Rayleigh — 
Benard convection arise even when the background density stratification is stable. 
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1 Introduction 



Double-component convection is relevant to as diverse fields as small-scale 
oceanography [1], astrophysics [2], geodynamo [3], crystal growth [4], colloidal 
suspensions [5], and soap films [6]. Convective flows are also commonly used 
for testing the ideas related to transition to turbulence and nonlinear pattern 
formation [7,8]. In addition, double-component flows where a distinction be- 
tween the components comes from component different boundary conditions 
are of basic significance for large-scale environmental phenomena. These phe- 
nomena range from Langmuir circulations [9] to the global ocean thermohaline 
circulation and climate change [10,11]. 
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One major aspect of double-component flows with different boundary condi- 
tions is that the effect of such conditions can be conceptually analogous to that 
of different diffusivities in conventional double-diffusive convection [12,13,14]. 
In particular, unequal diffusion gradients forming in perturbed state due to 
different boundary conditions trigger convection analogously to the classical 
double-diffusion. Such analogy has been introduced in [15,16,17,18] as a result 
of the generalization of an earlier idea highlighted in [19]. 

This work reports the existence of a novel mechanism of double-component 
convection. The nature of this mechanism stems from different boundary con- 
ditions but is not underlain by the differential diffusion caused by unequal per- 
turbation gradients of the components (differential gradient diffusion). This 
mechanism has been identified in a horizontal layer of pure, Boussinesq fluid 
where an unstable gradient of the component whose boundary values are fixed 
is combined with a stable stratification due to the flux boundary conditions for 
another component. It manifests itself in the form of flnite-amplitude steady 
flows before the onset of the respective linear instability. Not resulting from 
differential diffusion, such flows arise even when the net background strati- 
fication is stable. An appropriate perturbation could thus trigger convection 
in a broad range of parameters where such convection could not have been 
previously anticipated. 



2 The problem formulation and solution procedures 

Let the diffusivities of the components be equal and let a stable stratification 
due to the component with flux boundary conditions be combined with an 
unstable gradient of the component with the fixed-value conditions. The dif- 
fusivities are set equal to eliminate the effects of the classical double-diffusion, 
and thus to examine the effects of different boundary conditions separately. 
This is analogous to the approach used in previous studies of conventional 
double-diffusive convection. In most such studies, the components with un- 
equal diffusivities have not been distinguished from each other in terms of 
boundary conditions. Equal diffusivities could, besides, characterize two so- 
lutes, as claimed in [20]. The ratio of the viscosity to diffusivity (which is the 
Prandtl number, Pr) would then be different from that for temperature in 
water {Pr = 6.7) used in this work. This parameter, however, is not expected 
to have a qualitative effect on the novel phenomenon reported herein. Equal 
diffusivities could also be viewed as eddy coefficients [9,11], when the ratio 
between the viscosity and diffusivity is close to its value used below. 

The background gradients arc represented by the Rayleigh numbers Ra = 
ga\AT\(f/ and Ra^ = g(3\dS/dx\d'^/ = /iRa. Here, a; is the (dimensional) 
vertical coordinate, d is the width of the horizontal slot, AT is the (dimen- 
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sional) difference between the values of temperature (standing for the compo- 
nent with fixed- vahie boundary conditions) at the lower and upper boundaries, 
dS/dx is the boundaries-prescribed (dimensional) vertical derivative of solute 
concentration (standing for the component with flux boundary conditions), 
a is the coefficient of thermal expansion, /3 is the coefficient of the density 
variation due to the variation of solute concentration, g is the gravitational 
acceleration, u is the kinematic viscosity, and k = = is the difTusivity of 
both components. The bar means that the respective variable is dimensional. 

The configuration just described is illustrated in Fig. 1 as ^ = 0, 9{> in 
Fig. 1) is the angle between the direction opposite to the gravity and that 
of the across-slot coordinate axis. With = d"^ /dx^ + d"^ jdy^ ^ = 8/ Or, 
dx = d/dx, and dy = d/dy, the equations describing the two-dimensional 
problem in Fig. 1 in terms of streamfunction ip, t, and S are: 

Pr{dr + dxiidy - dyiidx)ALi> = 

(dxt - dxS)sine - (dyt - dyS)cose + PrA|z^, (1) 
Pr{d, + dx^dy - dyiPdx)^, = Al^„ 1 = 1,2. (2) 

Here .^i and ^2 stand for t and the across-slot velocity u = —dyip, the along- 
slot velocity v = dxip, and r is the time. This problem was studied for 6' = and 
the periodic boundary conditions with period A = A/d=2in the along- slot 
direction by continuation [21] of (finite-difference) steady solutions in Ra and 
H, the numerical approach was the same as in [15,17,18]. For clarification of one 
relevant issue arising for > (this case will be explicitly identified below), 
continuation in 9 was also used. S = was set at the middle points of the 
across-slot boundaries, along with the otherwise periodic conditions. Such a 
condition is needed to specify the scale of S and the phase of a nontrivial steady 
solution. The time evolution of a linear perturbation initially imposed on the 
steady flows was computed to examine the solution stability. The stability of 
the conduction base flow to steady disturbances was also analyzed for different 
wave numbers k = 27r/A. This was done by searching for the smallest- it!a 
singularity of the matrix resulting from the application of boundary conditions 
to the general solution of the steady, marginal linear stability problem. 



3 Background 

An inflnitesimal disturbance imposed on the conduction state (Fig. 1, 9 = 
0) would lead to the formation of unequal perturbation gradients. Due to 
the differential gradient diffusion in perturbed state, a rising (sinking) fluid 
element would experience the buoyancy force directed downwards (upwards) 
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[19]. The buoyancy force is thus expected to act against the sense of rotation 
of a small-amphtude perturbation celL This permits amphtude growth of the 
perturbation cells changing their sense of rotation with adequate frequency, 
as illustrated in [19] for the inviscid fluid. (As highlighted in [15,16], this 
effect makes the present configuration, as well as the configuration in [19], 
analogous to the diffusive regime of conventional double-diffusive convection 
[12,13].) Such oscillatory instability also arises in the viscous fluid (Fig. 2), 
when the stable flux stratification increases. A detailed discussion of the effect 
of viscosity on manifestation of the oscillatory instability on different scales is 
beyond the scope of this work. Its main idea is given in [22]. 

(In the framework of current discussion, an oscillatory perturbation could be 
viewed as a standing wave, i.e., as the convective cells whose sense of rotation 
changes periodically in time. The prescribed across-slot-boundary values of S 
can prevent traveling-wave disturbances from being detected in the present 
formulation. However, traveling waves are also expected to arise from such 
Hopf bifurcations as H in Fig. 2 if the translation symmetry of the conduction 
state is allowed for [23]. Their presence, in particular, may affect the stability 
of steady branches in Fig. 2.) 

Since the differential (gradient) diffusion in the perturbed state (Fig. 1, = 0) 
results in the growth of oscillatory infinitesimal disturbances, it is expected to 
oppose growth of the small stationary perturbations. (The sense of rotation 
of a stationary-perturbation cell does not change. The cell amplitude could 
thus grow only against the above effect of differential gradient diffusion.) In- 
deed, Ra'^{k) in Fig. 3 increasingly exceeds Ra^{k)/{1 — yu) as /i grows from 0. 
[i?a°(A;)/(l — ii) are the values Ra^{k) would take on if both boundary con- 
ditions were of the fixed- value type.] The deviation of Ra'^{k) in Fig. 3 from 
Ra^{k) / {1 — iJ.) is particularly pronounced for the large scales (small k), where 
diffusion is most effective. With respect to infinitesimal steady disturbances, 
therefore, fiux boundary conditions for the solute stratification stabilize the 
conduction state compared to the single- component problem with the unstable 
fixed- value gradient. 



4 Results and discussion 

The basic result of this work is that a novel physical mechanism due to different 
boundary conditions (Fig. 1, ^ = 0) gives rise to finite- amphtude convective 
steady fiows where the conduction state is stable to the infinitesimal steady 
disturbances (Fig. 2). As in convection resulting from differential gradient dif- 
fusion, one element of this mechanism is disparate responses of the component 
stratifications to convective motion. In the present mechanism, however, the 
feedback to convective perturbation arises from finite-amplitude Rayleigh — 
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Benard convection. This is essentially different from differential gradient diffu- 
sion in [15,16,17,18,19]. With such new feedback, different boundary conditions 
are found to result in a purely nonlinear manifestation of convection being due 
to the statically stable net vertical stratification. 

Finite-amplitude convective flows are illustrated in Fig. 4. As the convection 
amplitude increases [Fig. 4(a), (b)], the ratio of the across-slot solute concen- 
tration scale to such scale in the background state decreases. This is the result 
of an increasing number of solute isolines moving "outside" the flow domain, 
especially in the regions of across-slot motion. Such behavior is associated with 
the flux conditions permitting solute isolines to cross the boundaries. The re- 
spective ratio for the temperature, however, remains equal to one, even when 
convection becomes well-developed [Fig. 4(c), (d)]. The flxed- value boundary 
conditions maintain the vertical temperature scale by preventing the intersec- 
tion of isotherms with the boundaries. They also increase the thermal gradi- 
ent near the wall towards which the across-slot component of convection is 
directed. 

The unstable density gradients thus formed in the regions of across-slot mo- 
tion in Fig. 4 substantially exceed i?a°(7r) in Fig. 3(a). [i?a°(A;) in Fig. 3(a) 
represents the onset of Rayleigh — Benard convection.] Such gradients result in 
a horizontal density difference between two streamline points, as in Rayleigh — 
Benard convection. This gives rise to (positive) convective feedback that main- 
tains the disparity between component gradients. Such feedback forms even 
when the net background stratification is neutral or stable [Fig. 4(c), (d)]. In 
such cases of the present formulation, linear steady instability does not arise, 
as in the scenario first proposed in [24]. As /i increases, the cells at the smaller- 
amplitude branch, Al, change their form to utilize the regions with maximal 
gradient disparity more efficiently [Fig. 4(e)]. 

The analogy in the physics of oscillatory instability between the present con- 
figuration and the diffusive regime of conventional double-diffusive convection 
[13] does not seem to apply to flnite-amplitude steady instability in both 
these problems. The flnite-amplitude mechanism in [13] hinges on the dis- 
parity between component diffusivities. Such a disparity makes the unstable 
temperature gradient relatively insensitive to a convective perturbation. This 
gives rise to the feedback maintaining convection. In the present mechanism, 
differential gradient diffusion plays only a stabilizing role as the (feedback) 
unstable density stratiflcation arises from the interaction of perturbation with 
boundary conditions. In contrast to [13], in particular, subcritical steady con- 
vection arises in the present problem even when the density stratification is 
statically stable. 

Nonlinear Rayleigh — Benard convection also gives rise to feedback in the 
finite-amplitude steady instability in binary fiuid [25] , if the separation ratio 
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is negative. However, the binary-fluid finite-amplitude mechanism is underlain 
by the dependent nature of the stabilizing background (Soret) solute gradi- 
ent, rather than by boundary conditions. Such solute gradient is thus largely 
destroyed when its conduction-state relation to the unstable temperature gra- 
dient is relaxed by a finite convective perturbation. 

Let > in Fig. 1 and let ^ = 1. When = 7r/2, finite-amplitude steady fiows 
arise due to differential gradient diffusion [18]. For large enough convection 
amphtudes {Ra >~ 20000), the dissimilarity between the nature of such flows 
at 9 = 11/2 and that of finite-amplitude steady convection at ^ = was found 
to give rise to a region of hysteresis in 6. Within this region, different convective 
steady flows with A = 2 coexist. One of these solutions was continued from 6 = 
tt/2 in decreasing 9. For say Ra = 31000, it extends to 6* ~ 0.387r/2. Another 
such flow was continued from ^ in increasing 9. For Ra — 31000, it extends 
to ^ fti 0.5l7r/2. For the same Ra{^ 31000), such hysteresis phenomenon was 
not found between 9 = tt/2 and 9 = n, where steady convection is driven 
only by differential gradient diffusion [15,16,18]. A detailed analysis of the 
flow transformations will be reported separately. 



5 Conclusions 

Underlying a fundamentally new convective mechanism, different boundary 
conditions thus result in finite- amplitude steady convection arising in a hor- 
izontal layer of pure fluid well before the respective linear instability. Such 
nonlinear convective flows also exist when the vertical density stratiflcation 
is statically stable. The new mechanism has to remain basically relevant for 
Le 7^ 1, as well as under different Pr (as say for two solutes) and stress- free 
boundaries, among other changes. It also raises the issue of three-dimensional 
effects. Another issue is the existence of an analogous effect (as well as of the 
effects of differential gradient diffusion) when the buoyancy forces are replaced 
by the forces due to surface tension. All this leads to a new perspective on 
the role of convection and different boundary conditions in double-component 
fluid systems, including the large-scale systems relevant to global environmen- 
tal processes. 
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a-^ = 6S-6T 
t = -R a/2 Sx=-txRa 

t = Ra/2ySx=-|Li Ra 
No-slip walls 




Fig. 1. The problem in a horizontal {9 = 0) and inclined {9 > 0) slot. 6p = p — po 
is the variation of the (dimensionless) density, p, due to the variations SS and 
6T of solute concentration S and temperature T = (Ti + T2)/2 + t with respect 
to their reference values, at which the density is po; Ti and T2 are the boundary 
temperatures, a = gd^/nu. Pr = v/k = 6.7, Le = ht/hs = 1; and ks (= n) are 
the component diffusivities. The fluid is of the Boussinesq type. 
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|i=0.6 |i=0.4 [1=0.2 \i=0.l 




10000 Ra 20000 



^|i=1.2 



Fig. 2. 9 = 0. Schematic structures of the steady flows with minimal period A = 2 
for /i G [0, 1.2] (Pr = 6.7, Le = 1). $ is an abstract measure of the steady flows 
that distinguishes between different steady solutions, specifies the location of the 
singularities (limit points and bifurcation points), and represents the flows arising 
from a symmetry-breaking bifurcation as a single branch. The background states 
are depicted by the horizontal lines with arrows (for fi = 1.2, this is the coordinate 
axis). The solid lines stand for the stable solutions. The dashed lines represent the 
flows being unstable to either steady or both steady and oscillatory disturbances. 
The dotted lines stand for the solutions being unstable to oscillatory disturbances 
alone. [Instability of the background state to steady disturbances with wave number 
= 27r (A = 1), arising before that with = tt (A = 2) at = 0.6, is not shown.] B 
is the bifurcation standing for the steady linear stability boundary for wave number 
A; = TT (A = 2); it changes its criticality at jj. just below 0.2 and moves to infinite 
Ra as fi ^ 1. L is the limit point standing for the finite-amplitude steady stability 
boundary for the flows with minimal period A = 2. It moves to Ra « 25772 as 
increases to 1.5. Al and A2 are the unstable and stable (to steady disturbances) 
branches associated with the limit point, respectively. H is a Hopf bifurcation. For 
/i = 1, it arises at Ra 29000. 
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Fig. 3. 9 = 0. Curves of marginal linear stability to steady disturbances for different 
// {Le = 1), i?a^(fc). (a) = 0; (b) = 0.1; (c) = 0.2; (d) /x = 0.3; (e) fi = 0.4; 
(f) fi = 0.5; (g) fi = 0.6; (h) /x = 0.7; (i) /x = 0.8; (j) n = 0.9. 
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Fig. 4. 9 = 0. Convective steady flows (Pr = 6.7, Le = 1); X = 2. streamlines; 
S: isolines of solute concentration; T: isotherms {t in Fig. 1). The solid and dashed 
streamlines represent the clockwise and counterclockwise rotation, respectively. The 
actual values of >S' and t are equal to 10^ times the respective values in the figure, 
(a) 12 = 0.6, Ra = 12046, branch Al (directly unstable); (b) = 0.6, Ra = 6500, 
branch A2; (c) /j. = 1, Ra = 24000, branch A2; (d) /x = 1.2, Ra = 24000, branch 
A2; (e) /x = 1.2, Ra = 50090, branch Al (directly unstable). 
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